arXiv:1505.01222vl [cs.IT] 5 May 2015 


1 


Joint Interference Alignment and Bi-Directional 
Scheduling for MIMO Two-Way Multi-Link 

Networks 

A. M. Fouladgar, O. Simeone , O. Sahin , P. Popovski and S. Shamai (Shitz) 


Abstract —By means of the emerging technique of dynamic 
Time Division Duplex (TDD), the switching point between uplink 
and downlink transmissions can be optimized across a multi-cell 
system in order to reduce the impact of inter-cell interference. 
It has been recently recognized that optimizing also the order 
in which uplink and downlink transmissions, or more generally 
the two directions of a two-way link, are scheduled can lead 
to significant benefits in terms of interference reduction. In this 
work, the optimization of bi-directional scheduling is investigated 
in conjunction with the design of linear precoding and equaliza¬ 
tion for a general multi-link MIMO two-way system. A simple 
algorithm is proposed that performs the joint optimization of the 
ordering of the transmissions in the two directions of the two-way 
links and of the linear transceivers, with the aim of minimizing 
the interference leakage power. Numerical results demonstrate 
the effectiveness of the proposed strategy. 

Index Terms - MIMO, Two-way communications, Scheduling, 
Interference alignment, Linear precoding. 


I. Introduction 

Two-way communication is one of the most common modes 
of operation for wireless links, particularly for cellular and 
Device-to-Device (D2D) systems (see Fig. []]). The conven¬ 
tional approach to the design of a system comprising multiple 
interfering two-way links consists of two separate phases: 
at first, one fixes a transmission direction independently for 
each link; and, then, the physical-layer parameters, such as 
powers or beamforming vectors, are optimized in a centralized 
or decentralized fashion so as to maximize some system- 
wide performance criterion. This is, for instance, the standard 
approach for cellular systems, in which the scheduling of 
uplink and/or downlink transmissions is performed on a per- 
cell basis as a preliminary step. 

The emerging technique of dynamic Time Division Duplex 
(TDD) breaks with the conventional approach discussed above 
in that the switching points between uplink and downlink 
transmissions in a frame are optimized jointly across all cells 
based on the current channel conditions 12-0. The key 
motivation for this paradigm change is the observation that 
the selection of the duration of the uplink and downlink 
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Figure 1. A network of two-way MIMO communication links. Examples 
include uplink-downlink cellular links and two-way D2D channels. 


transmission phases in a cell can have a significant impact 
on the interference observed by the interfered cells. 

With dynamic TDD, the ordering of the transmissions in the 
two directions of a two-way link is fixed. However, the inter¬ 
ference configuration depends significantly on such scheduling 
decisions. An illustration of bi-directional scheduling is pro¬ 
vided in Fig. [2] for two links: in each frame comprising two 
slots of fixed duration, each link can operate in the left (L) 
to right ( R ) direction, i.e., L —f R, in the first slot and then 
in the opposite direction, R —> L, in the second slot, or it 
can schedule first the R -A L direction and then the L —t R 
direction. As in we refer to this binary choice for each 
link as the interference spin of the link. 

Reference m studies the optimization of bi-directional 
scheduling in the presence of fixed switching times for single¬ 
antenna links. In this work, instead, we investigate the in¬ 
terplay of multi-antenna transmission and reception with bi¬ 
directional scheduling for fixed switching times and a general 
MIMO two-way multi-link system. This study is motivated by 
the observation that the choice of the interference spins across 
the network defines the channel matrices of the interfering 
channels on which the linear transceivers operate. Therefore, 
the capability of multiple antennas to mitigate interference 
depends on the scheduling decisions in a fundamental way. 
Algorithms are proposed that perform either the separate or 
the joint optimization of the interference spins and of the 
linear transceivers, with the aim of minimizing the interference 
leakage power and hence obtaining enhanced interference 
alignment solutions in the sense of (7). 

Notation: We use lower case fonts for scalars and uppercase 
bold fonts for vectors and matrices. I,/ represents the d x d 
identity matrix and 0 mxn is used to indicate the M x TV 
zero matrix, tr (A) denotes the trace of the matrix A and 
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Figure 2. Illustration of the definition of interference spin Si for each link 
i. A frame consists of two consecutive slots as shown by bold lines. 


A h is the conjugate transpose of matrix A. ^ in (-) returns 
a truncated unitary matrix spanning the space associated with 
the d smallest eigenvalues of the argument matrix. <C MxN 
represents the set of complex-valued M x N matrices and 
IJA/xAr represents the set of truncated unitary matrices with N 
orthonormal columns. X is the complement of binary variable 
X. The indicator function 1[-] returns one if the argument is 
true and it returns zero otherwise. 

II. System Model 

We consider a wireless network consisting of K MIMO 
two-way interfering links. An illustration is shown in Fig. Q] 
for a scenario with two cellular links and a D2D link. Each 
two-way link fc is equipped with Ni,k antennas at the left- 
hand node Lk and with N Rjk antennas at the right-hand node 
Rk . Note that the labeling of one end of each link as “right” 
or “left” is arbitrary. Each link operates using Time Division 
Duplex (TDD), and two-way communication takes a frame 
consisting of two successive slots: in the first slot, each fcth 
link may operate either in the direction Lk —t R k , so that the 
Lk node is the transmitter and the Rk node is the receiver, or, 
vice versa, in the direction Rk —> Lk, and the successive slot 
is used in the opposite direction. All links are synchronous 
and are assumed to be always backlogged so that there is a 
continuous stream of frames. 

Based on the discussion above, each link k operates in 
either direction Lk —> Rk or Rk —> Lk in even slots and in 
the opposite direction in odd slots. Following we define 
the order in which the two directions are scheduled as the 
interference spin or, for short, spin of a link. Specifically, the 
fc-th link is said to have a 0-spin if it operates in the direction 
Lk —>• Rk in the odd slots and the direction Rk —> Lk in 
the even slots; otherwise, when it operates in the direction 
Rk —»• Lk in the odd slots and in the direction Lk —» Rk 
in the even slots, the fc-th link is said to have a 1-spin. We 
refer to Fig. [2] for an illustration. The spin of the /;:-th link 
is denoted by Sk £ {0,1}. In order to simplify the notation, 
in the following, we set Nj Jt k = N R and N Rt k = Nr for all 
k £ 1C, although the generalization is straightforward. 

A. Signal Model 

Referring to Fig.[3]for an illustration of the main definitions, 
the signal y R k £ received by the fcth receiver when 

active in the Lk —> Rk direction is given by 

y Rk = h R k ,L k *L k + ^2 (i [Sj = s k ] u RkyLj tc Lj (i) 

jeic\{k} 

+1 [Sj ^ Sk] U RkyRj x R .) + z Rk , 
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Figure 3. Illustration of the definition of the MIMO channels. 


where the matrices U RkRj £ i[^n r xn l anc j H Rk;Rj £ 
C NrxNr re p resen t the channel responses from the node L ;j to 
Rk, and from the node Rj to Rk, respectively; the vectors x/ 
and x. R . are the signals transmitted by the nodes Lj and Rj, 
respectively; and the vector z Rk is the additive noise, which 
is distributed as z Rk ~ CAf(0, a 2 1) and is independent across 
the link index fc. The first term in © is the signal received 
from the desired transmitter while the second term represents 
the interference from all the other links that have the same 
spin as link fc and the third accounts for the interference 
caused by links that have opposite spin. Similarly, the signal 
y Rk £ C Nlx1 received by the fcth receiver when active in the 
Rk —> Lk direction can be written as 

y L k = n Lk , Rk X Rk + (! i S 3 = Sk] Hl^AjXRj (2) 

je>c\{k} 

+l[Sj^S k ]U Lk ^ Lj )+z Lk , 


with analogous definitions. 

The channel matrix ~iA Rk j., and similarly for the other pairs 
of nodes can be written as H Rk , Lj = cr RktLj il Rk<L:j , where 
the parameter a Rkt L j models path-loss and shadowing over 
the link Rk — Lj in both directions R k —»• Lj and Lj —> Rk 
and is defined as 


a Rk , Lj 



v 


(3) 


where D re f is a reference distance, D RkRj is the distance 
between node Rk and node Lj, rj is the path loss exponent 
and ^ Rk ,L } is the log-normal shadowing component, which 
is distributed as t; Rk ,Lj ~ W(0, z/ 2 I) and is independent for 
different pairs of nodes; and H RktR . is an N R x N R channel 
matrix that accounts for the effect of fast fading on the channel 
Lj —» Rk and is assumed to have i.i.d. CAf(0, 1) entries. The 
channel matrices may or may not be reciprocal, so that, in 
general we have the inequality H # fcikj f- H|^. Rk . 

The channel matrices are assumed to be constant for at 
least a frame. The analysis below will be limited to any 
period including frames in which the channel matrices remain 
approximately constant. 

When transmitting, nodes fc/. and R k , first perform chan¬ 
nel coding, producing codeword vectors s Rk £ C dk x1 and 
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s R k £ <C dkXl , which are distributed as S£ fc ~ CJ\T(0,Pk/dkT) 
and s ji k ~ CAf(0 : Pk/dkl), respectively, where parameter 
dk < min (Nr, N^) represents the number of data streams 
exchanged on the two-way link and Pf~ is the received average 
signal-to-noise ratios, (SNRs), excluding shadowing effects, at 
the reference distance D K f. Note that we assume that the power 
and the number of data streams is the same in both directions 
for each link, although generalizations are straightforward. 
We also assume for simplicity that equal power allocation 
is performed across all data streams. After channel coding, 
node Lk performs linear precoding with a truncated unitary 
precoding matrix £ \] NLXdk yielding the transmitted 

baseband signal 

x Lfc =V L)i s Lk (4) 

and similarly node Rk performs linear precoding with a 
truncated unitary precoding matrix VR k £ \] N n xdk yielding 
the transmitted baseband signal 


Note that the defined interference covariance matrices depend 
on the interference spin variables Sk, k £ /C, and on the linear 
precoding and equalization matrices. 

III. Bi-Directional Scheduling and Interference 
Alignment Optimization 

In this section, we propose algorithms that perform the 
optimization of the interference spins vector S = [Si,..., Sk] 
and of the linear transceivers U = [{Ufl fe }^L 1 , {Ui fc }^L 1 j 
and V = [{V fl jf =1 ,{V L jf =1 ]. Following 0, in order to 
approximate the interference alignment conditions, we adopt 
as the optimization criterion the interference power leakage. 
Note that 0 only considers one-way communication links 
and hence the optimization therein is limited to the linear 
transceivers. Specifically, we formulate the problem of interest 
as the minimization of the total received interference leakage 
power as 


x R k = VR k SR k 


(5) 


On receiving the baseband signal y n k , node Rk performs 
linear equalization with a truncated unitary equalization matrix 
£ U NRXdk yielding the output signal 

so = URjft fc! (6) 

and similarly, on receiving the baseband signal y L k , node Lk 
performs linear equalization with a truncated unitary equaliza¬ 
tion matrix U L k £ \] NLXdk yielding the output signal 


(7) 


SR k = V L k yL k - 

The covariance matrix of the interference signals that affects 
the equalized signal s n k , is given as 


qRx _ ttH 
'“iR k ~ U i?j 




i — dj 


'i = Sk] , Lj V Lj V )L . 

+1 fa ± S k ] (8) 

while for the equalized signal s L k we have 

El 

— dj 
£/C\{fc} 3 


Q'o = u"( 


+1 [Sj ± S k ] H 
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V, V ' TT 


H TT H 
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(9) 


It is also useful to define the covariance matrix of the overall 
interference that is caused by node Rk when transmitting as 


q«; = V; 


">( £§( 1[s ' 


+1 [Sj ± S k ] H 


i = Sk] H" ^U l ,U"H 


H 




V 


Rtf 


Lj j Rk 


( 10 ) 


while the covariance matrix of the overall interference that is 
caused by node Lk when transmitting is given as 


Q™ = <1 E ^[Sj=Sk]n^ Lk v Rj v%u 

V'e/c\{fc } 3 


+i [Sj ± s k ] nl tLk v Lj v^.u Lj , Lk ))v Lk . 


H ■ 


Rj }L k 


( 11 ) 


minimize 

s,u,v 


E 


tRX . t RX 
fcG/C 1 R k ' 1 L k ’ 


( 12 ) 


or, equivalently, as 

mmiimze E fceK El +EEr- ^ 13 ) 

where we have the implicit constraints that S £ (0,1} A hold 
and that the matrices U and V be truncated unitary with the 
mentioned dimensions; moreover, the quantities 

^ = tr(Q**) and/^ = tr(Q^) (14) 

measure the received interference power for link k in the two 
directions, while 


% = ^(QS) and/-=tr(Q-) (15) 

measure the interference power caused by the transmitters of 
link k in the two directions. 

We observe that it is possible to extend the proposed 
approach to other performance criteria including the max- 
SINR method of |7j], but this will not be further pursued here. 
Instead, the focus is on approximating the interference align¬ 
ment conditions as per the interference leakage minimization 
(ILM) scheme in 0. We start by discussing an algorithm 
that performs the separate optimization of the interference 
spins and of the linear transceivers in Sec. IIII-AI and then we 
introduce an algorithm that carries out the joint optimization 
of interference spins and linear transceivers in Sec. IIII-BI 


A. Separate Optimization 

A first simple solution is that of first determining the spin 
variables S and then perform linear transceiver optimization 
using the ILM algorithm as in 0. 

Spin Optimization: In order to optimize the spin variables, 
we propose to solve the problem of (IT2l >. or equivalently (IT3l) . 
by setting all precoding and equalization matrices equal to 
the identity matrix in (ISll- dUl) . This optimization can either 
be carried out by exhaustive search in the space (0, 1} R or 
by using the algorithm proposed in 0. Moreover, in order 
to minimize the overhead associated to the selection of the 
vector S, one can perform the optimization of S based only on 
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long-term CSI, namely path loss and log-normal shadowing. 
This enables the spins to be updated only at the time scale of 
the long-term fading variability. In this case, problem (IT2l) . or 
m, is tackled by setting all matrices U and V to the identity 
matrix and by averaging the interference powers over the fast 
fading variables. This leads to the problem 


min J2 


kGlC 


lE 


1 [Sj=S k \ 


\jeic\{k} 

+1 [Sj^Sk] 

1 [Sj=S k ] 


D 


n f \ v ^Rk' L i 

ref ' 10—ro ~N l N r 


Rk,L. 


Dre f 


D 


Rk,Rj 


• 10 ^ Nfi 


n f x 7 * 

ref ' IQ^^NlNr 


+l[Sj?S k ] 


D L k ,R 

Aef ^ 


D Lk 


10 15 Nr 


(16) 


which can be tackled as explained above. 

Interference Leakage Minimization (ILM): After find¬ 
ing the optimized spin vector S = [Si, S 2 , Sk] as ex¬ 
plained above, the linear transceiver matrices V and U can 
be calculated using the ILM method of J7] Alg. 1] using 
full CSI. The ILM algorithm aims at minimizing the in¬ 
terference power over a A -user MIMO interfering channel 
with generic channel matrices {Gfcj € C NkXM:i } fc over 

the precoding matrices {P^. € U MfcXdfc }^_ and decoding 
matrices {Dj, eO NkXdk } k _ 1 as illustrated in Fig. [4] The 
algorithm runs over a given number Ailm of iterations and 
is based on the calculation, which can be implemented in 
a decentralized fashion, of the leading eigenvalues of given 
covariance matrices 0. We denote the outcome of the ILM 
algorithm as 




=ILM 


A,r 




(17) 


where we do not make explicit the dependence on the initial 
conditions for simplicity of notation. We recall the ILM 
algorithm for convenience in Table Algorithm |T] 

Based on the given vector S, we apply ILM separately to 
the even and odd slots. For the even slots, we assign the 
corresponding channel responses matrices H RktR ,, H Rk R , 
H L k ,R } and H Lk>Lj to G k j such that 

G kd = l[S i =0]l[S j = S k ]H Rh , Lj 
+ 1 [Sj = 0] 1 [Sj / Sk] H Lk , Lj 
+ l[S j = l]l[S j = S k ]U Lh , Ri 
+ l[S J = l]l[S J ^S k ]U RkMj (18) 


and 

P fe = 1 [Sk = o] V Lk + 1 [Sk = 1] V Rk (19) 


and 

D k = 1 [S k = 1] U Lk + 1 [S k = 0] U Rk (20) 

for all j. k £ 1C. Instead, for the odd time slots, the assignment 
above are obtained by substituting S with its entry-wise 
complement S. 


Algorithm 1 Interference Leakage Minimization |j7) 

1: Start with initial precoding matrices P y with P^P, = I 
for all j £ K: . 

2: Begin iteration with setting i = 1 

3: Compute interference covariance matrix at the receivers k £ 
1C: 

QP - E (L,PP' , G' , ; 

j&IC\{k} 

4: Compute the linear equalization matrix at each receiver k £ 

D* = (Qf) 

5: Compute interference covariance matrix for each transmitter 
j with j £ 1C: 

Qf = E G k,P^k'Gk.j 

keic\{j} 

6: Compute the linear precoding matrix for each transmitter j 
with j £ 1C'. 

Pj = 4„ (Qf) 

7: If * > Ailm, exit; otherwise i = i + 1 and go back to 3. 



Transmitter K Receiver K 

Figure 4. Illustration of a A'-user MIMO interfering channel. 


B. Joint Optimization 

In this section, we propose a technique that carries out the 
joint optimization of interference spins and linear transceivers 
based on full CSI. In order to avoid the exponential complexity 
of exhaustive search, which scales as 2 K , the proposed scheme 
explores a subset of possible spin vectors guided by a simple 
local search criterion. Specifically, at each step i, the algorithm 
explores a new spin vector S^) based on the previous vector 
S(* _1 ) by flipping the spins in S^ -1 ) corresponding to the 
minimum number of links, causing the largest transmitted 
interference powers. More precisely, we define the set > 
to include all spin vectors explored up to iteration (i— 1) and as 
Iop t the minimum interference leakage power observed up 
to iteration (i — 1) along with the corresponding spin vector 
SopE- Moreover, we define as /f - 1 11 and /f ^ lf the 
interference powers ( fl5l ) calculated using ILM as discussed 
in the previous subsection for the current spin vector S^ -1 ). 
At each step i, we flip the spin variables of the vector 
S(* _1 ) in the order of decreasing caused interference power 
if + /f until a new spin vector ^ 
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Algorithm 2 Joint Spin and Linear Transceiver Optimization 
1: Initialize the vector S^and set = 0, 


and IIp to a very large value. 

2: Begin step i = 1. 

3: Run two instances 

ILM ({G fclj }J . =1 , {dj}f =1 , N 1LM , {P,}f =1 , {D,}f =1 

with (IT8] >-(l2()t. one instance with S^ -1 ^ in lieu of S and 
one with S''~ l ;i in lieu of S. Calculate the resulting 
total two-way transmitted interference leakage, i.e.. 


of 


r TX (i-1) _ tTX (i-1) 
1 k ~ 1 R k 

leakage, i.e., I^ x 
k G /C. 


rTX (i-1) 


and received interference 


= I' 


RX (i-1) tRX (i-1) 

+ h 


Rk 


for all links 


4' If V 7 TX 1 -r-, ^ 

11 2^k 1 k T ± k — -‘opt 


t RX (i-1) T (i- 1) 


set /opt = /(* ^ and 


Sopi = S (l 1} ; otherwise, set I™ = 7 0 ( * t 11 and S^{ = S, 


opt 


x opt 


J opt 


a(i-l) 

J opt 

5: Update the set of explored spin vectors U 

6: Flip the smallest number of spin variables of the vector 
S (l_1 ) in the order of decreasing caused interference power 


7™ ^ + lj x 11 until a new spin vector ^ S <1 ^ is 


obtained. 

7: If i > F, stop, otherwise set i = i 


1 and go back to 4. 


is obtained. The variables S^, 7™ lJ X 7^ and Sop t 
are updated accordingly based on ILM for the vector S^. We 
perform a number F of comparison steps with F < 2 K . 

Remark 1. At each step i, rather than using the largest 
transmitted interference powers criterion, we could gener¬ 
ate a new vector S h> by flipping the spin of the link, or 
links, with the largest received interference powers, namely 
1 R + l Lk , when using the previous spin vector 

S (' — 1 1 . Extensive numerical results, discussed in part in the 
next section, reveal that the transmitted interference power 
criterion is generally to be preferred. 

IV. Numerical Results 

In this section, we evaluate the performance of the consid¬ 
ered techniques via numerical results. We consider a 100m x 
100m area in where all the links are located randomly. We 
set rj 3 and £ = 5dB for all links. There are K\ 
shorter-distance links with a transmitter-receiver distance of 
Di = 25 m and 7\2 longer-distance links with distance of 
D 2 = 50m, so that K = K\ -f K 2 - The shorter-distance 
links may represent D2D connections and the longer-distance 
links may represent WLAN channels between a device and 
an access point. In keeping with this interpretation, it is also 
assumed that both nodes of a shorter-distance link have the 
same average SNRs at the reference distance D re f = I)\ of 
Pl = 1’i-t = lOdB, while for the longer-distance links we 
have Pl = 30dB and Pr = 50dB. Note that the corresponding 
average SNRs at the distance D 2 are /’/, UFd'/ P 2 ) n = 21dB 
and Pn(D re f/D 2 ) v = 41dB. The number of antennas and the 
number of data streams are the same for all the nodes, namely 
TVr = Nn = 4 and d = 2, respectively. The number of ILM 
iterations is IVilm = 20. 



Figure 5. Average relative interference leakage power (dB) versus the number 
of comparison steps for the proposed approach and for reference schemes 
= 3, — 2, Nl = Nr = 4, d = 2, and -/Vhjvi = 20). 


Fig. 0 illustrates the average interference leakage power 
versus the number of comparison steps of the proposed Al¬ 
gorithm 0 with K\ = 3 and K 2 = 2. For comparison, we 
show the performance of exhaustive search, of ILM run on a 
randomly selected spin vector and of the separate optimization 
schemes discussed in Sec. IIII-AI For Algorithm 0 the spin 
initialization is done using the result of separate optimization 
using long-term CSI, i.e., by solving problem The 

average interference leakage is normalized to the average 
interference leakage power obtained with random spins. The 
performance of the proposed techniques is seen to be very 
close to that of exhaustive search as long as around 10 com¬ 
parison steps are performed. This is substantially smaller than 
the required number of comparison steps in exhaustive search, 
i.e., 2 5 = 32. This advantage is less pronounced if one uses 
instead the received interference power criterion, as discussed 
in Remark 0 Moreover, the interference leakage power can be 
significantly reduced with respect to the considered baseline 
techniques, e.g., by more than lOdB over a random spin 
selection. This demonstrates that the interference alignment 
capabilities of the system are enhanced by the adoption of joint 
optimization. Finally, we observe that separate optimization 
based on long-term CSI perform very similarly to separate 
optimization based on full CSI. 

Fig. [6] shows the average interference leakage versus the 
number of shorter-distance links with K 2 = 2 and the number 
of comparison steps set to 10 for the proposed algorithm. 
The interference leakage power is normalized to the average 
interference leakage power obtained with random spins. While 
the interference leakage power increases with K\, the relative 
performance of all schemes remains fairly constant, except for 
the algorithm based on the received interference power (see 
Remark 0, which, as seen above, tends to require a larger 
number of comparison steps when Ki is large enough. 

Fig.0shows the average sum-rate, which is calculated using 
standard capacity formulas (see, e.g., Q), versus the number 
of links 7\ i with K 2 = 0. The relevant gains of optimizing the 
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Figure 6. Average relative interference leakage power (dB) versus the number 
of shorter-distance links K\ ( K 2 = 2, Nl = Nr = 4, d = 2, and TVilm — 
20 ). 



Figure 7. Average sum-rate versus the number of links K\ (K 2 = 0, 
Nr = Nr = 4, d = 2, and A/jlm — 20). 


also seen that the interference leakage gains do not necessarily 
translate into sum-rate gains, for which separate optimization 
with long-term CSI appear to be sufficient to reap most of the 
performance advantage. Interesting future work includes the 
investigation of a distributed version of the proposed scheme. 
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spin vector are confirmed to hold also in terms of achievable 
two-way rate. Separate optimization is, however, shown to reap 
most of this advantage, even when using only long-term CSI. 
Finally, we note that increasing the number of links K\ is at 
first beneficial as long as interference is manageable, but it 
leads to a decrease of the rate as soon as the number of links 
is large enough (here, larger than 9). 

V. Concluding Remarks 

In this paper, we have investigated a wireless network with 
MIMO two-way interfering links that operate using TDD 
with fixed switching times. We have proposed to optimize 
the bi-directional schedule, i.e., the order of the transmission 
directions, along with the linear transceivers, with the aim of 
minimizing the interference leakage power. Numerical results 
show that the performance of the proposed technique in terms 
of interference leakage is very close to that of a scheme based 
on exhaustive search, but at a substantially lower complexity, 
and that significant gains can be obtained as compared to a 
random selection of the interference spins. Interestingly, it is 





